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We study interacting Rashba-Dresselhaus fermions in two spatial dimensions. First, we present a 
new exact solution to the two-particle pairing problem of spin-orbit-coupled fermions for arbitrary 
Rashba and Dresselhaus spin-orbit interactions. An exact molecular wave function and the Green 
function are explicitly derived along with the binding energy and the spectrum of the molecular 
state. In the second part, we consider a thermal Boltzmann gas of fermionic molecules and compute 
the time-of- flight velocity and spin distributions for a single fermion in the gas. We show that the 
pairing signatures can be observed already in the first-moment expectation values, such as time-of- 
flight density and spin profiles. 
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I. INTRODUCTION 

The concept of bound states of fermion pairs arises 
ubiquitously in physics. The interior of a neutron star is 
believed to be a superfluid of bound neutron pairs held 
together by an attractive component of the inter-nucleon 
interactions. 1 In particle physics, mesons are a bound 
state of two fermions, a quark and an anti-quark, with 
the attraction provided by the strong interaction. Ar- 
guably the most well-known two-fermion bound state in 
condensed matter physics is the Cooper pair 2 , bound to- 
gether in some cases through phonon exchange 3 , which 
is the key component in the understanding of supercon- 
ductivity. 

A bound state of a fermion pair also plays an im- 
portant role in the field of ultracold gases of neutral 
Fermi atoms, where atoms in two different hyperfine 
states can be used to realize pairing. With the advent 
of the Feshbach resonance 5 9 experimentalists are now 
able to adjust both the magnitude and sign of the two- 
body scattering length by simply tuning an external mag- 
netic field 10 14 . This has opened doors to a new field 
which explores the crossover between a BCS superfluid 
of weakly bound Cooper pairs and a Bose-Einstein con- 
densate (BEC) of tightly bound fermionic molecules. 15 18 
The BEC of such molecules have been observed in time- 
of-flight (TOF) experiments which show a characteris- 
tic bimodal distribution similar to the one observed in a 
BEC of bosonic atoms. 19 23 

In addition to the unprecedented control over the 
strength and nature of inter-atomic interactions exper- 
imentalists are also capable of fine-tuning with relative 
ease the structures and parameters of these atomic sys- 
tems using quantum optical techniques. These devel- 
opments open the possibility of emulating various solid 
state systems using ultra-cold atomic systems, and us- 
ing them to gain insights into the outstanding problems 
in strongly correlated condensed matter physics. 4 ' 24 ' 25 In 



this regard, the theoretical proposals 26 31 for generat- 
ing synthetic Abelian and non-Abelian gauge fields for 
a gas of neutral fermionic atoms, and the recent ex- 
perimental realization of the phenomenon by the NIST 
group 32 34 , brings great excitement to the community. 
While a spatially inhomogeneous Abelian gauge field can 
mimic the role of an external magnetic field, a uniform 
non-Abelian gauge field is capable of introducing spin- 
orbit (SO) interactions akin to Rashba and Dresselhaus 
couplings found in semiconductor electronic systems. 30 
These achievements may serve as a key ingredient when 
realizing some of the most intriguing phases of matter in 
condensed matter physics, including fractional quantum 
Hall liquids 35 ' 36 and topological insulators 37 ' 38 . 

A fundamental feature of low-density electronic sys- 
tems with Rashba SO interaction is the enhanced ten- 
dency for bound pair formation due to a SO-induced in- 
crease in the density of states of low-energy electrons. 39 
An inspection of the single-particle energy spectrum for 
a Rashba particle reveals a ring of minima in energy in 
contrast to a single minimum at the momentum-space 
origin for a particle without SO coupling. Indeed, it 
has been noted, for instance, that SO interaction can 
be remarkably beneficial for superconducting pairing. 40 
The experimental realization of SO interactions in ultra- 
cold Fermi gases has motivated many recent theoretical 
works. 41 Some of these include the BCS-BEC crossover 
physics 42 46 and the effects of population imbalance and 
superfluidity 47 52 in the presence of SO interactions, as 
well as properties of the BEC ground state in the strong 
gauge-field limit 53 . 

In this work, we study tightly-bound molecules of 
Rashba-Dresselhaus fermions in two spatial dimensions, 
and investigate the properties of a low-density thermal 
gas of such molecules. We first develop in detail the ex- 
act solution to the problem of two Rashba-Dresselhaus 
fermions interacting via an attractive short-ranged s- 
wave interaction. The wave function for the bound 
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molecular state is then obtained by solving for the bound 
state energy and the molecular energy spectrum. We 
then focus on a low-density Boltzmann gas of these bound 
molecules and compute the single- fermion density ma- 
trix for one of the two fermions that form a molecule 
in the gas. With this density matrix we evaluate the 
fermion velocity and spin density distributions in mo- 
mentum space that can be inferred from a TOF experi- 
ment. We compare and contrast these distributions with 
those corresponding to a non-interacting Boltzmann gas 
of SO-coupled fermions. The main finding that stems 
from this consideration is that the pairing of a fermion 
with its molecular partner has an experimentally observ- 
able imprint on these distributions, and that these im- 
prints are absent when fermions are not interacting. We 
show that these signatures of pairing appear already, for 
instance, in the first-moment density expectation value 
(n(k)). Higher-order correlations as a means to probe 
many-body states has previously been proposed in the 
context of ultra-cold atoms. 54 

The paper is organized as follows. The exact solution 
to the above-mentioned interacting two-particle problem 
is presented in Sec. II. The properties of the bound 
molecular state are explored in Sec. Ill: analytic ex- 
pressions for the binding energy and the molecular spec- 
trum in different physical limits are shown in Sees. Ill A 
and IIIB, respectively. Effects due to a Zeeman term is 
briefly discussed in Sec. IIIC. In Sec. IV, we consider 
a thermal gas of non-interacting molecules and focus on 
the single-fermion density matrix for such a system (Sec. 
IV B). There, we compute the TOF velocity (Sec. IV C) 
and spin (IV D) distributions for a single fermion in the 
system. In Sec. IV E we discuss how these distributions 
can be measured experimentally. 



II. THE EXACT TWO-PARTICLE WAVE 
FUNCTION 

Let us begin with a system of two "spin" -1/2 fermions 
with two-dimensional Rashba and Dresselhaus SO inter- 
actions moving in two spatial dimensions. The single- 
particle Hamiltonian is given by 

H RD (k) = e k I+a (& x k y - a y k x )+f3 (a x k x - cr y k y ) , (1) 

where ek = k 2 /2m with m being the mass of the par- 
ticle, & = (a x ,a y ) are the usual Pauli matrices and 
k = (k xi k y ) is the particle momentum. We will be using 
units where h = 1 throughout. The term proportional 
to a (/?) is known as the Rashba (Dresselhaus) SO in- 
teraction. The hats denote 2x2 matrices acting on the 
"spin" -space of the particle. Here, the "spin" degree of 
freedom may be a synthetic degree of freedom. Hamilto- 
nian (1), for instance, can be realized in ultra-cold atomic 
systems where atoms with multiple internal levels move 
in the presence of a spatially modulated laser field. 30 ' 63 
For the two-particle system the kinetic energy contri- 



bution to the Hamiltonian is a 4 x 4 matrix given by 

#o = iWki) ® I + I ® H RD {k 2 ), (2) 

where k^ with i G {1,2} is the momentum of the i-th par- 
ticle and the check indicates an operator acting on the 
two-particle Hilbert space. For some interaction Hamil- 
tonian i?i n t the total Hamiltonian of the system is then 
given by H = H + # in t- 

The energy eigenstates of our two-particle system are 
four-component spinors and are labeled by two indices, Q 
and n, where Q = ki + k2 is the centre-of-mass momen- 
tum and n labels both the various bound and scattering 
states. Any energy eigenstate |-0 n (Q)) can be written as 
a superposition of momentum eigenstates as 

\MQ)) = J ^xn(Q,k)|Q,k), (3) 

where k = (ki — k2)/2 is the relative momentum. 
Note here that Xn(Q?k) are four-component spinors 
owing to the different possible spin states of the two- 
particle system. We will be using the triplet-single ba- 
sis, \S,s z ) = {|1,+1),|1,-1),|10),|00)}, to express the 
spinor (see (20)), unless otherwise stated. The stationary 
Schrodinger equation is then 

[ff-IB n (Q)]|^ n (Q)) = 0, (4) 

where E n (Q) is the energy of the (Q,n) state. Inserting 
(3) into (4), we obtain 

[ff (Q,k)-iE n (Q)] Xn(Q,k) 

+ J (k|ff int |k / )xn(Q,k , ) = 0. (5) 

We note that a general isotropic interaction in two di- 
mensions can be written as 

oo 

<k|ij int |k')= Y, vHMV (0k - 0k,) , (6) 

l= — oo 

where 0k is the angle between k and the x-axis. We will 
hereafter assume that the interaction is short-ranged. At 
low energies and long wavelengths, scattering amplitude 
is dominated by the contribution of the s-wave scattering 
as long as the relative momentum k satisfies kR e <C 1, 
where R e is the characteristic radius of the interaction. 
Due to the anti-symmetry of the two-particle wave func- 
tion we must project out the singlet component of the 
wave function for the s- wave component (i.e. 1 = 0). 
Dropping higher order harmonic contributions, we may 
then write the potential as (k| Hi nt |k') — >> Vq^ s \ where 
^( s ) = 1 00) (00 1 is the singlet projector. We are then led 
to rewriting (5) as 

[tf (Q,k)-l£ n (Q)] Xn(Q,k) 

— ^) x „(Q,k') = 0. (7) 
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Let us now define the inverse Green function as 
^(Q.k) = H (Q,k) - LE„(Q) and the function c% 



100) = J |^« Xn (Q,k'). 



Then we obtain the expression for the spinor, 
X„(Q,k) = -c^y C(Q,k)|00). 



(8) 



(9) 



The task now is to invert the inverse two-particle Green 
function, Sf n " 1 (Q, k). Hereafter, we will refrain from writ- 
ing the index n explicitly since it enters only to label the 
eigenenergies. We now introduce the single-particle in- 
verse Green functions for particle i, namely, g~ 1 (k i , Q) = 
H RD (ki) -i(£7 n (Q)/2) = iy(k<, Q) +acr • (b< x z), where 
y(ki,Q) =e k , -E n (Q)/2 and 



b i; = 



1 7 
7 1 



(10) 



Here, j = /3/a measures the relative strength between 
the Rashba and Dresselhaus couplings. The two-particle 
inverse Green function is then 

^- 1 (Q,k)= [r 1 (ki,Q)®I + i®r 1 (k 2 ,Q)] • (11) 

g _1 (k^, Q) can be diagonalized using the unitary matrix 

iii = exp j-i^ (m • 6") | , (12) 

where = b^/ft^. We then obtain ^j# _1 (k^, Q)i^ = 
Ii/ n (ki, Q) + a^a 2 . Therefore, £f _1 (Q, k) can be inverted 
using the composite unitary transformation U = Ui (£) U2 
and we obtain 



S?(Q,k) = ^(Q,^ 



(13) 



with 
£(Q,k) 

and 



Diag{[D 1 (Q,k)]- 1 ,[D 2 (Q,k)]- 1 , 

[^(Q.krMAiQ.k)]- 1 }, (14) 



£>i(Q,k) = s + a(6i + 6 2 ) 
D 2 (Q,k)=s + a(6 1 -6 2 ) 
D 3 (Q,k) = s-a(6 1 -6 2 ) 
D 4 (Q,k) = s-a(6i+6 2 ). 



Here, we have defined 



s(Q,k) 



E n (Q), 



(15) 
(16) 
(17) 
(18) 



(19) 



An explicit expression for the Green function matrix (13) 
is presented in Appendix A. Inserting (13) into (9) the 



spinor %(Q, k) can be expressed in the triplet-singlet ba- 
sis as 



x(Q,k) 



cqVq 
"d(Q,k) 



/ <(Q,k) 
-*?*(Q,k) 
<(Q,k) 
V ^ s (Q,k) 



(20) 



The denominator d(Q, k) is given by 



d(Q,k) =s 4 -4a 2 s 2 



(21) 



where B = bi +b2, b = (bi — b2)/2. The coefficients for 
the triplet and singlet components are explicitly given by 

*tf(Q, k) = zV^ [5a6e- i0b - a 3 (B • b)Be _i * B ] (22) 
(Q, k) = 2isa 2 [B x b] z (23) 
^ s (Q,k) = s(s 2 -a 2 £ 2 ), (24) 

where 0& = tan _1 (6 2y /6 x ) and similarly for (/>#. Note 
that both «e^(Q,k) and «e^(Q,k) are odd functions of 
the relative momentum k, while <£/ s (Q,k) is even. This 
ensures that the wave function %(Q, k) is anti-symmetric 
under interchange of two particles: Q — >> Q, k — >> — k 
and <Ji o CT2- (20) together with the Green function 
(Al) represent the exact solution, and is the main new 
technical result of this work. 

The normalization constant cq can be obtained from 
the orthonormality condition for the energy eigenstates 
|V(Q)>, i.e. (V(Q)|V(Q)} = (QIQ) for all Q (and n), 

|cq|2= [ v ° 2 J (0(oo|^(Q,kMQ,k)|oo)] \ 

(25) 

The solution is complete once the energy spectrum 
E n (Q) is obtained for all the eigenstates. 



III. PROPERTIES OF THE BOUND 
MOLECULAR STATE 

If an energy eigenstate (Q,n) describes a bound state 
of our two fermion system its spectrum must satisfy 
E n (Q) < 2E m i n for some values of Q. Here, E m i n is 
the minimum value in the single-fermion spectrum for 
Hamiltonian (1), which is 



2m 



±afe v / l + 7 2 + 27sin(20 /c ), (26) 



where <pk = tan _1 (/c 2/ //c x ). We find £? m in = —ma 2 (l + 
7) 2 /2. 7 = /3/a, again, is the relative strength between 
the Rashba and Dresselhaus couplings. We may now 
define the momentum-dependent binding energy A n (Q) 
through E n (Q) = 2E min (l + A n (Q)/ma 2 ). The condi- 
tion for a bound state then translates to A n (Q) > for 
some Q. 
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A. The binding energy at Q = 

Owing to the oddness of the triplet coefficients in (20) 
with respect to k, one finds that only the singlet com- 
ponent survives once one integrates over k in (8). This 
then leads to the eigenvalue equation for the bound state 



d 2 k 

(2^) 2 



X (Q,k)=c Q |00). 



(27) 



Making use of (9) leads to the self-consistency condition 
1 r cPk 



V 



(2n)< 



(00|<f(Q,k)|00>, 



(28) 



which, together with (1), (2) and (Al), has an explicit 
form 



/ 



d 2 k s(s 2 -a 2 B 2 ) _ _J_ 
(2tt)2 d(Q,k) " V ' 



(29) 



Introducing the dimensionless variables, 
Q k A 



ma 



ma 



5 n 



ma* 



v = mV , (30) 



and the dimensionless energy variable £ = k 2 the self- 
consistency condition for q = reduces to 



d 2 ^ 



£-e n (0) 



v Jo (2^) 2 (£-e n (0)) 2 -4a(^- 



(31) 



Here, we have defined the dimensionless energy spectrum 
en(q) = et(l + 5 n (q)), where e t = 2e min = -(1 + j) 2 is 
the dimensionless threshold energy and a(<p) = 1 + j 2 + 
27 sin 2(j). We have also introduced the UV cutoff A ~ 
1 / (maR e ) 2 which is set by the characteristic radius of the 
interaction potential R e . The typical 5 n (0)-dependence 
of the right hand side of (31) is plotted in Fig. 1 for 
7 = 0, 7 = 0.1 and 7 = 0.5. We find that for any given 
attractive interaction vo < 0, there is a single bound 
state. We label this state by n = 0, but refrain from 
explicitly writing the index. For the isotropic case (i.e. 
7 = 0) and for < 5(0) <C 1, we find that the binding 
energy is given by 



8(0) 



NT 
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=: 5 R (0). 



(32) 



In the weakly anisotropic regime (i.e. < 7/5(0) <C 1) 
we find the lowest 7-corrections to be 



5(0) « 5 R (0) 



7 



'6 R (0) 



7 



5 R (0) 



(33) 



while in the strongly anisotropic regime (i.e. 7/5(0) ^> 1) 
we obtain 



8(0) 



A 



(l + 7) : 



_ 87r^7 

-e hoia+7). 



(34) 




8(0) 



FIG. 1: Plot of the right hand side of (31) as a function of 
the binding energy 5(0) > 0. Here, we have chosen 7 = for 
the solid line, 7 = 0.1 for the dashed line and 7 = 0.5 for the 
dot-dashed line. While an algebraic divergence is observed 
as 5(0) — > in the 7 = case, a logarithmic divergence is 
observed in the presence of the Dresselhaus coupling (7 > 0) . 
The UV cutoff is A = 100. 



Here, we have assumed A ^> 1 + 7. We remind the reader 
that 7 = /3/a is the relative strength between the Rashba 
and Dresselhaus couplings. Similar calculations for the 
binding energy have been done previously for various SO 
coupled Fermi gases. 42 46 However, the consideration of 
arbitrary two-dimensional Rashba and Dresselhaus cou- 
pling strengths in two spatial dimensions and the pre- 
sentation of the exact wave function, to the best of our 
knowledge, has not been done. 

The crossover in the binding energy from an algebraic 
dependence on the interaction strength (c.f. (32)) to ex- 
ponential (c.f. (34)) is directly related to the crossover 
in effective dimensionality for bound state formation. 39 
Let us say that a SO coupled system in D spatial dimen- 
sions has a set of single-particle minimum-energy states 
in momentum space that forms a d-dimensional surface. 
The effective dimensionality for bound state formation 
is given by D e R = D — d. As can be seen from the 
single-particle spectrum (26) the isotropic case possesses 
a ID manifold of minimum states (d = 1) while once 
the Dresselhaus coupling is finite, this manifold of mini- 
mum states is reduced to two points in momentum space 
(d = 0). The effective dimensionalities for the two cases 
are then D e R = 1 and D e ff = 2, respectively. The ex- 
ponentially small binding energy (34) and the effective 
dimension D e ^ = 2 corresponding to that case is consis- 
tent with a bound state problem in quantum mechanics 
of a particle moving in a potential well in two dimensions. 

The binding energy can also be computed in the regime 
where 5(0) ^> 1. Since 7 ~ 1, we will simply consider 
7 = 0. If the cutoff scale is still the largest scale, such 
that 5(0) <C A, we find the binding energy to be 



5(0) ^Ae i«oi. 



(35) 
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with the corresponding effective mass of 



o 

q y -o.05 



FIG. 2: Numerical results for the dimensionless molecular 
spectrum e(q) = E(Q)/ma 2 as a function of the dimension- 
less molecular momentum q. (a) 7 = and |^o| -1 — 2 are 
used, (b) 7 = 0.1 and |^o| _1 = 1.5 are used. Anistropy in 
the molecular effective mass can be observed in (b) for finite 
Dresselhaus coupling (c.f. (38)). The plot shows the bound 
state ceases to exist for momenta larger than the critical mo- 
mentum q c rsj y/S(0). 



B. The molecular spectrum 

We first provide an approximate analytic expression 
for the spectrum valid for q 2 /S(0) <C 1, from which the 
effective mass of the molecule can be extracted. In the 
weakly isotropic limit, the binding energy to &(q 2 ) reads 



(q) « 6(0) 



8 



7 



2 S R (0) 



sin(20 g ) 



M eff (7) 



4m(l + 7) 2 



1 + 67 + 7 2 + (1 - 7) 2 sin(2<^) ' 



(39) 



(36) 



Interestingly, we find that isotropy in the effective mass 
is restored for equal Rashba and Dresselhaus interac- 
tion strengths (7 = 1), the so-called persistent spin 
helix point 55 57 . Let us now consider the effective 
mass of a free fermion, m e ff(7), near the bottom of 
the single-particle spectrum (26). We find that the 
molecular effective mass and the single- fermion effec- 
tive mass are related in a one-to-one relation, namely, 
M e ^(j) = 2m e ff(7). The dimensionless molecular spec- 
trum e(q) = E(Q)/ma 2 can be straightforwardly ob- 
tained through the relation e(q) = e t (l + 5(q)). 

The molecular spectrum can be obtained numerically 
for arbitrary values of q. The results for the spectrum as 
a function of q are plotted in Fig. 2 for 7 = and 7 = 0.1, 
respectively. \vo\ was chosen to be 0(1) such that we are 
in the regime of 6(0) <C 1. We see the anisotropy in 
the spectrum for 7 = 0.1 (Fig. 2(b)). We also see that 
for large enough momenta, q > q c , a bound state ceases 
to exist. The scale for the critical momentum is set by 
q c ~ ^/S(0). The existence of this critical momentum 
tells us that the fermions do not remain bound once the 
molecular kinetic energy exceeds the binding energy. 

An estimate of the binding energy can be made if a 
gas of these fermionic molecules is rotated. We consider 
the gas confined to the x?/-plane and a rotation about 
the z-axis. For simplicity, we will take the pure Rashba 
case here, i.e. 7 = 0. In the dilute limit, where the gas 
can be treated as a classical (Maxwell-Boltzmann) gas, 
the velocity field under the rotation ft = Qz is given 
by a rigid rotation v(R) — QR<fi, where <fi is the unit 
tangent vector on the xy-pldine and R = \J x 2 + y 2 is the 
radial distance from the axis of rotation. The critical 
momentum scale then introduces a critical distance scale 
R C1 where, for R < R C1 the molecules are still bound but, 
for R > R C1 we have a gas of unbound fermions. If we 
introduce a dimensionless distance r = RQ/a the critical 
distance is given by r c = g c /4 ~ s/5(0). Therefore, the 
radial distance at which the two phases of bound and 
unbound fermions meet gives an estimate of the binding 
energy. 



where <j) q = tan -1 (q y / f q x ) . For 7 = 0, the isotropy in 
the spectrum is restored as expected. The corresponding 
effective mass for the molecule is then given by 



M eff ( 7 ) ^4m 



(37) 



In the strongly anisotropic regime, the binding energy to 
ff(q 2 ) reads 



o(q) « 6(0) 



q 2 1 + 67 + 7 2 + (1 - 7) 2 sin(20,) 
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(1 + 7) 4 



. (38) 



C. Effects due to a Zeeman field 

We briefly discuss the effects of a synthetic Zeeman 
field, H = (H x , H y , H z ), which couples to the pseudo- 
spin of the fermions. We introduce a Zeeman term, 
Hz = —(H x a x + HyCFy + H z cf z ) to the single-particle 
Hamiltonian (1). In Appendix B the self-consistency con- 
dition is rederived for a Zeeman field oriented in a general 
direction in R 3 . Here, we will explicitly consider the field 
pointing in the z-direction, and for the pure Rashba case 
where 7 = 0. We find that the self-consistency condition 
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at q = is given by 



A. Qualitative discussion of the result 



v 



Jo 2tt 



(£ - e (0)) 2 - 4h 2 



2ir (^- eo (0))[(e-e (0)) 2 -4e-4/i2]' 

(40) 

where the dimensionless Zeeman field h = H/ma 2 . Since 
the threshold energy for molecular formation is now e t = 
— (l+/i 2 ), the dimensionless energy for the bound state at 
q = reads e (0) = -(1 + h 2 )(l + (5(0)). For h<h c = l, 
(40) gives 



1-h 2 



1 



vTTT? 



(41) 



where again the approximation holds for 6(0) <C 1 and 
A ^> 1. For h < /i c , a bound state exists for all |^o|, 
but the binding energy approaches zero as h —> h c . Once 
h > h c , the integral on the right hand side of (40) be- 
comes bounded as 6(0) — >■ and, therefore, a bound state 
ceases to exist for sufficiently weak attractive interaction 
strengths. This introduces a quantum phase transition 
at a critical coupling Vq separating phases with bound 
and unbound fermions. 



IV. LOW-DENSITY RASHBA-DRESSELHAUS 
MOLECULAR GAS 

Let us now consider a low-density gas of TV tightly- 
bound Rashba-Dresselhaus molecules confined to two 
spatial dimensions. We assume that the gas is equili- 
brated at some temperature T, which satisfies Tbkt <C 
T«A. Here, Tbkt is the temperature at which the gas 
undergoes a Berezinskii-Kosterlitz-Thouless phase tran- 
sition into a superfluid, and A is the molecular binding 
energy. Presumably, in this temperature range, the gas 
is in the dilute limit, nA T <C 1, where n is the areal 
density of the molecules and At = h/ ' yJlnJtkBT ', with 
jtft := min{M e ff(7)}, is the mean thermal wavelength 
evaluated for the smallest molecular effective mass. The 
gas may then be modeled as a thermal Boltzmann gas of 
uncondensed molecules. 

A molecule can interact with other atoms and 
molecules in the gas. Indeed, scattering between atoms 
and bound molecules as well as between two molecules 
was considered in depth in many works. 58 62 Here, we as- 
sume that the gas is sufficiently dilute so that we may ne- 
glect atom-molecule and molecule-molecule interactions 
to first order. 

Our aim first is to obtain TOF velocity and spin distri- 
butions for a single fermion atom for the molecular gas. 
These distributions can be inferred from a spin-resolved 
TOF experiment, where both the trap and the SO cou- 
pling are turned off and the depaired fermions are allowed 
to expand freely. We then contrast these distributions to 
corresponding distributions for a Boltzmann gas of un- 
bound (non-interacting) fermionic atoms and discuss the 
striking differences between the two cases. 



For the gas of non- interacting molecules ("gas A") 
single-fermion distributions must be extracted from the 
single-molecule density matrix. For the gas of non- 
interacting fermion atoms ("gas B") the correspond- 
ing distributions are obtained directly from the single- 
fermion density matrix. For gas A the single-molecule 
density matrix is given by 



QA 



Tr{e-P**y 



(42) 



where the Hamiltonian H = Hq + H int was given in (2) 
and (6) in Sec. II. For gas B, the single-fermion density 
matrix reads 



-PH RD 



Qb 



(43) 



where the Hamiltonian was given in (1). The key differ- 
ence in the single-fermion distributions for gases A and B 
stems from the fact that while the single-fermion momen- 
tum eigenstate for gas B, |k), is an eigenstate of Hrd, 
the two-fermion momentum eigenstate, |Q,k), for gas A 
is not an eigenstate of H. For gas B the momentum op- 
erator commutes with the Hamiltonian, and the velocity 
distribution for the Boltzmann gas is trivially given by 



F(k) = (k| g B |k)oce-^, 



(44) 



where was given in (26). For gas A, we first obtain 
the diagonal elements of the molecular density matrix, 
(Q, k| qa |Q, k), and the velocity distribution is extracted 
by tracing out one of the fermions. Since |Q, k) is not an 
eigenstate of H, we must introduce the energy eigenstates 
|VVi(Q)) m order to replace H by its expectation value. 
This results in the diagonal elements which schematically 
has the form 

F(k 1 ,k 2 ) = (Q,k|^|Q,k) 

KX , (Q,k) X (Q 1 k)e-W, (45) 

where E(Q) is the molecular spectrum obtained in Sec. 
IIIB and x(Q? k) is the spinor wave function obtained in 
(20). The extra factor involving x m (45) indicates that 
each fermion is correlated with its partner fermion due 
to interactions. In the free fermion case (c.f. (44)), such 
correlations are clearly absent. 

The correlations between fermions in gas A has an im- 
print on the various momentum distributions and make 
them distinct from the corresponding distributions in gas 
B. Moreover, these distinctions can be made within 
first-moment expectation values (e.g. (n(k))). These 
differences can, in principle, be inferred from TOF ex- 
periments. For small binding energies (i.e. 6(0) <C 1) 
the qualitative shapes of the various momentum distribu- 
tions, in general, appear similar for both gases. The key 
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difference arises in the width of the peak features found 
in the distributions. For free fermions (gas B) these peak 
features should have a gaussian profile with the width set 
by temperature, as it is clear from (44). In contrast, the 
peak features for gas A has a square-Lorentzian profile 
with the width set by the binding energy scale. In fact, 
for gas A, the peak features are determined by the wave 
function %(Q,k), and not the Boltzmann factor, which 
primarily determines the centre-of-mass distribution of 
the molecules (c.f. (45)). 

For r := k B T/ma 2 <C 5(0), the width in the peak 
structures for gas A should be markedly broader than 
those for gas B. The peak widths for gas B can be ad- 
justed by changing the temperature while such variation 
should not occur for gas A. On the other hand, increas- 
ing the attractive interaction strength, thus increasing 
the binding energy, while holding the temperature fixed 
should lead to a broadening in the width of the peaks for 
gas A only. 



B. Quantitative results for various distributions 

We now provide detailed calculations of the single- 
fermion distributions for gas A described above. Recall 
that the (dimensionless) molecular dispersion was given 
by e(q) = e t (l + 5(q)), with q-dependent binding ener- 
gies (36) or (38) which are valid for q 2 <C S(0). The valid- 
ity of (36) or (38) is ensured in our temperature regime, 
since r = kBT/ma 2 <C 5(0), such that essentially all of 
the fermions are bound and occupy single-molecule states 
very close to q = 0. The diagonal elements of the single- 
molecule density matrix (42) can then be approximately 
written as (see Appendix C) 

|c q | 2 (00 1 Syt( q , K )g( q , K ) | QQ) e -e(g)/r 
Jq,* K\ 2 <00| Sft(q, K)S?(q,K) |00> e~<^ ' 

where J q K = J j^j^- Here again, q and n are 
dimensionless momenta and q = K\ + k<i and k = 
(k\ — K2)/2. The velocity distribution for a single fermion 
can be obtained by integrating out the other, 

P(*l) = J ^P(«l,«2). (47) 

Note that P(k,\,k,2) = P(k>2,k>i) so one could have in- 
tegrated out either one or the other electron and would 
have arrived at the same probability. 

The single-fermion spin distributions can also be ob- 
tained analogously. The three spin operators for a 
fermion are Si = <3^/2, where i = x,y,z. The i-th com- 




FIG. 3: Single-fermion velocity distribution for the molecular 
gas for three different relative SO interaction strengths: (a) 
7 = 0; (b) 7 = 0.005; and (c) 7 = 0.5. Binding energy 
5(0) = 0.05 is used, and the temperature is taken to be small 
such that r < 5(0). 

ponent of the spin density is then given by 

:= (Si) (#61) = J ^| 

|c q | 2 <00| #(q, K)[Sj ® t]Sf (q, k) |00) e-<^ 
/ qK |c q | 2 (00|^t(q,«Mq,«)|00)e-^)/- ' 

Note here that the matrix Si <g> I acts in spin space 
spanned by the basis \cr u a 2 } = {|tt) , , lit) > 
We find that ^(^1) = 0, which is expected from solv- 
ing a single-fermion quantum mechanical problem with 
Hamiltonian (1). 

At low temperatures but still well above Tbkt, the 
Boltzmann factor in (47) and (48) is a strongly peaked 
function at q = with contributions becoming exponen- 
tially small for q > y/r. In the low temperature regime 
we are considering (r <C 5(0)) the remaining factors in 
(47) and (48) give a slowly-varying function of q on the 
scale of y/r. We may then drop the q-dependence in those 
factors and replace it by 0. As a result, the integral over 
a^2 can be done trivially and we arrive at 

P(k) w C (00| #(0, «)Sf (0, k) |00) , (49) 
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FIG. 4: Single-fermion spin densities in momentum space for 
7 = 0. S x (k) and S y (n) are plotted in (a) and (b), respec- 
tively. Binding energy 5(0) — 0.05 is used, and the tempera- 
ture is taken to be small such that r <C £(0). 



FIG. 5: Single-fermion spin densities in momentum space for 
7 = 0.5. S x (k) and S y (n) are plotted in (a) and (b), respec- 
tively. Binding energy 5(0) = 0.05 is used, and the tempera- 
ture is taken to be small such that r <C 5(0). 



where C is a normalization constant. Similarly, the spin 
densities can be approximately written as 



Si(K) « c (oo| ^ T (o, K)[Si iqsf (o, k) |oo) . (50) 



C. Velocity distribution 

Plots of P(pv) are shown in Fig. 3 for three different 
values of 7: (a) 7 = 0; (b) 7 = 0.005; and (c) 7 = 
0.5. The q = binding energy was taken to be 5(0) = 
0.05 and r = 10 -6 . We find that the maxima in the 
velocity distributions occur for values of momenta where 
the minima in the single-particle spectrum (26) occur. 
For 7 = the spectrum has a ring of degenerate minima 
at k = 1, and this is reflected in the distribution in Fig. 
3(a). The Dresselhaus interaction breaks this degeneracy 
and the spectrum yields two minima at (K x ,K y ) = ±(1 + 
7, 1 + 7)/y / 2- The two peaks in Fig. 3(c) coincide again 
with the locations of these minima. The distribution in 
the crossover regime between the isotropic and strongly 
anisotropic limits is plotted in Fig. 3(b). 

The fact that large weights are observed at spectrum 
minima is expected as we are in the dilute and low tem- 
perature limits where most of the fermions are occupying 
momentum states near the band minima. Although cor- 
responding distributions for free fermions (gas B) will be 
qualitatively similar to the distributions obtained here 
for gas A, they are quantitatively different. For the pure 
Rashba case (7 = 0), for instance, the width of the ring 
of maxima for gas B is set by the temperature (c.f. (44)). 



In particular, the peak profile has a gaussian profile with 
the width set by y/r. In contrast, for gas A, a radial cut 
of the velocity distribution is essentially proportional to 
X f (°^)x(0^) and yields 



P(n) oc 



(^ 2 + l + (5(0)) 2 +4^ 2 



[(^ + l) 2 +J(0)] 2 [(/^-l) 2 +J(0)] 2 

(2 + £(0)) 2 +4 
[4 + 5(0)] 2 [(^-l) 2 + J(0)] 2 



(51) 



for 5(0) <C 1. Therefore, the ring of maxima has a square- 
Lorentzian profile with the width set by y/5(0). In the 
limit 5(0) ^> r, we would thus expect the peak width to 
be much broader for gas A than the corresponding dis- 
tribution for gas B. If 5(0) is gradually increased while 
keeping the temperature fixed the ring of maxima for 
the molecular gas will get progressively broader. Once 
5(0) ^> 1, the ring can no longer be resolved and one 
obtains a single broad peak centred at k — 0. This is 
illustrated in Fig. 6(a) where we have evaluated (49) for 
5(0) = 5 (c.f. (35)). This evolution of the velocity dis- 
tribution as 5(0) is increased is unique to that molecular 
gas, for, in the case of free fermions, the ring of maxima 
should remain sharp with the width set by y/r. 



D. Spin density distributions 

Spin densities in momentum space, S x (n) and S v (k), 
are plotted in Figs. 4 and 5 for the regime r <C S(0) <C 1. 
We consider 7 = and 7 = 0.5 in Figs. 4 and 5, re- 
spectively. In analogy to the velocity distribution, the 




FIG. 6: Single-fermion velocity and S x distributions with 7 = 
in the large binding energy limit (i.e. 5(0) ^> 1); here, we 
have taken 5(0) = 5, and the temperature is still taken to be 
small such that r < 1. 

spin-densities obtained here appear very similar to the 
corresponding spin-densities obtains by solving a single- 
fermion quantum mechanical problem with Hamiltonian 
(1). However, again as in the case of the velocity distri- 
bution, the widths of the maxima are set by the bind- 
ing energy scale. In Fig. 6(b), S x density for the pure 
Rashba case is plotted for 5(0) = 5. There, the half- 
rings of maxima and minima, as seen in Fig. 4(a), are no 
longer resolved due to the large binding energy. 

E. Detecting the distributions in TOF experiments 

These momentum distributions will be directly observ- 
able through a spin-resolved TOF measurement. 32 The 
TOF signature will be dependent on the SO scheme used. 
In what follows we assume the effective SO coupling is 
induced using the TV-level scheme, with N = 4 63 , which 
can be implemented in the alkalis such as 6 Li. In this 
scheme, the dressed states have the form 

\Da) = \jZ^ aj/2 \^) ( 52 ) 

where |j) = e* Kj * r \j) is a bare hyperfine state, 
boosted by Kj = ma[— smfrj /2)e x + cos(7rj/2)e 2y ]. The 
two pseudo-spin states are given by |t) = \D\) and 
II) = \D*). 

The velocity-spin distributions shown in Fig. 4 and 5 
can be inferred through a spin-resolved TOF measure- 
ment. Such a measurement will give the velocity distri- 
bution = (I £ti e --(— 'W/2 \ Da ) (Ai ,|). At 
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1 Hk. 




FIG. 7: Single-fermion velocity distribution for the four bare 
spin states |1) , |2) , |3) , |4) are plotted in (a), (b), (c) and (d), 
respectively, for 7 = in the small binding energy regime (i.e. 
5(0) <C 1). Here, 5(0) = 0.05 is used, and the temperature is 
taken to be small such that r <C 5(0). 

low temperatures the state \Ds) and {D^) will not be pop- 
ulated, and the expectation value reduces to = 

Kit) (tl + II) (II + e i7rj/2 It) (II + e" i7rj/2 II) (II). These 
states have the equivalent momentum distributions 



<P1>(") = 






(53) 


<P2>(«) = 




- 2S x (k)) 


(54) 






+ 2Sy(K)) 


(55) 


{V 4 ){k) = 




+ 2S x (k)) 


(56) 



where Vj = \j) (j\ is a projective measurement into the 
bare spin state \j). It is therefore possible to recon- 
struct the P(k), S x (k) and S y (n) momentum distribu- 
tions from a TOF measurement. 

These spin distributions are plotted in Figs. 7 and 8. 
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(a) 




FIG. 8: Single- fermion velocity distribution for the four bare 
spin states |1) , \2) , |3) , |4) are plotted in (a), (b), (c) and 
(d), respectively, for 7 = 0.5 in the small binding energy 
regime (i.e. 6(0) <C 1). Here, 6(0) = 0.05 is used, and the 
temperature is taken to be small such that r <C 6(0). 

In the pure Rashba limit, the momentum distribution of 
the four bare states have the same structure up to a ir/2 
rotation in momentum space. In the anisotropic limit, 
7^0, the bare spins have the same two-peak struc- 
ture, but with different relative amplitudes. Similar to 



the momentum distribution, the bare spin distribution 
can be distinguished from the case of non-interacting SO 
coupled fermions by the dependence of the width of the 
distribution on the interaction strength. 



V. CONCLUSION 

We investigate the properties of a low-density molecu- 
lar gas of Rashba-Dresselhaus fermions in two spatial di- 
mensions. The gas is considered at sufficiently high tem- 
peratures such that it can be considered as a thermal gas 
of uncondensed tightly bound molecules. The description 
of the gas is based on the exact solution to a quantum 
mechanical problem of two Rashba-Dresselhaus fermions 
interacting via an attractive short-ranged s-wave interac- 
tion. We compute the single-fermion density matrix for 
the gas and evaluate the fermion velocity and spin distri- 
butions. By making comparisons to corresponding dis- 
tributions for a Boltzmann gas of free fermions we show 
that these various distributions can be used to probe pair- 
ing of fermions in the molecular gas. Moreover, we find 
that the signatures of pairing appear in first-moment ex- 
pectation values. We discuss a spin-resolved TOF mea- 
surement from which the various distributions can be 
inferred in an experiment. This result is not specific 
to SO coupled systems considered in this work. Anal- 
ogous signatures of correlations should appear generally 
in single-fermion distributions for a system composed of 
interacting fermions. 
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Appendix A: Explicit expression for the Green 
function matrix (13) 

Recall that the Green function matrix was given by 
£f(Q,k) = UD(Q, k) U\ where U = u\ 0^2, and ui and 
i)(Q,k) were defined in (12) and (14), respectively. We 
then find 



/ (+ + ++) 



Sf(Q,k) = 



i e -i4>2 ( + _ + _) ie -i4>i (+ + ) 



1 + 02 



_i e .*n( + _ + _) 



(++++) 



-ie v 



i{<j>i+<p2) 



L(_|__| ) e *(01-02)( + 1_) 



'(+--+) *e-** 1 (+ + — ) 



++) 



j e -#2(_|_ _ _|__) 



(Al) 



(+--+) -ie** 1 (+ + — ) -ie 1 ^ (+-+-) 



++) / 



where fa = tan 1 (b Vj i/b Xj i), and 



Here, pi = ±1 and Di were defined in (15)-(18). 



(P1P2P3P4) = Pi Di ' 



(A2) 
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Appendix B: Self-consistency condition with a 
Zeeman field 

In this appendix, we give a brief derivation of the self- 
consistency condition in the presence of a Zeeman term 
of the form Hz — H <r. where a = (a,a z ). The inverse 
single-particle Green function in this case can be written 
as g _1 (lq, Q) = fi?/(k^ Q) + a& • (hi x z) — H z cr zi where 



where 
Bi = 

e 2 = 



1/ 


IT + a 2 b 


2 -W 


H - a 


2 b 


2 -H 2 






H'-ab 





H'+ab H'-ab 



(B6) 
(B7) 



bi = bi + -H' =: (b^,6 z i ), 
a 



(BI) 



where H' = (—H Xl —H yi H z ). It can be diagonalized 
with a unitary transformation 



exp 



7r + 2s in- 1 ^]K- ( T)|, (B2) 



where n- = b-/6-. By the composite unitary transfor- 
mation U = u\ <g> U2 the two-particle inverse Green func- 
tion can be written as £f(Q,k) = UD(Q, k)£A, where 
Z)(Q,k) is identical to (14) but with ^ replaced by b^. 
The two-particle wave function coefficients can be formed 
as before (c.f. (9)). The self-consistency condition then 
reads 



1 



S 



d 2 k 



2a 2 (j? + &y fl _ h 



d(Q,k) 

(B3) 

where B = bi + b 2 and b = (bi — b 2 )/2, and the de- 
nominator is defined as before (c.f. (21)) but with b and 
B replaced by b and B, respectively. f\ and / 2 are given 

by 



h = 



a 



2hb 2 
-26i6 2 ©2, 



h 1 - 



+m 



©i 



(B4) 
(B5) 



Appendix C: Obtaining P(ki,k,2) 

P(ki, K2) is given by the diagonal elements of the den- 
sity matrix (42). Since the momentum eigenstates |q, k) 
are not eigenstates of the Hamiltonian, we must insert 
a resolution of unity using energy eigenstates. We then 
obtain 

(q,K|e-^|q,K> = 

£ / (q,K|Vn(q'))(V'n(q , )|q,«)e- e " (q ' )/T . (CI) 

At low temperatures, we can safely assume that most of 
the fermions are in the n = bound state. Therefore, we 
may approximate the above matrix element by 



(q,«|e-^ ff |q,« 



f (q ) K|Vo(q'))(V'o(q , )|q,«)e- eo(q ' )/T . (C2) 
Using (3) and (9) we may write (C2) as 



(q,K|e-^|q,K)^|c q |V 2 (q|q) 

x (00| % f (q. ^)%(q, k) |00) e - eo(q)/r . (C3) 

(46) follows directly from this expression. 
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